CUBATURE ON WIENER SPACE: PATHWISE CONVERGENCE 



CHRISTIAN BAYER AND PETER K. FRIZ 

Abstract. Cubature on Wiener space [Lyons, T; Victoir, N.; Proc. R. Soc. Lond. A 8 
January 2004 vol. 460 no. 2041 169-198] provides a powerful alternative to Monte Carlo 
simulation for the integration of certain functionals on Wiener space. More specifically, 
and in the language of mathematical finance, cubature allows for fast computation of Eu- 
ropean option prices in generic diffusion models. 

We give a random walk interpretation of cubature and similar (e.g. the Ninomiya- 
Victoir) weak approximation schemes. By using rough path analysis, we are able to estab- 
<-h lish weak convergence for general path-dependent option prices. 
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1. Introduction 



Cubature on Wiener space (Kusuoka [6,7], Lyons and Victoir [9], see also Litterer and 
Ph Lyons [8 ], Ninomiya and Victoir |[T3l . Ninomiya and Ninomiya [ 12|) provides a powerful 

Ph . alternative to Monte Carlo simulation for the integration of certain functionals on Wiener 

space. As of present, these functionals are of the form f{Sr) where St is the image of a 
(i-dimensional Brownian motion under the Ito-map (the solution map to a stochastic differ- 
ential equation); the aim of cubature on Wiener space is then to provide a fast numerical 
algorithm to compute E [/(St)], where the expectation is taken over the ^/-dimensional 
Wiener measure. 

In the language of mathematical finance, cubature deals with European option prices in 
generic diffusion models. Although some exotic options can be handled in this framework 
(e.g. Asian options, by enhancing the state-space) general path-dependent options are not 
included in the presently available analysis on cubature methods. It must be admitted that 
t^J" . cubature has been designed for fast evaluation of payoffs of the type f (St)', but even so, it 

_±I ' may maintain its benefits in mildly path-dependent situation and, in any case, convergence 

(^) , to the correct value will be considered a minimal requirement by most users. 

CO ' The answer to "How can it fail to converge to the correct value?" is not trivial: cu- 

bature methods are essentially derived from replacing Wiener-measure but a path-space 
measure supported on smooth paths {oi,-}, subject to certain technical conditions relating 
to the iterated integrals of these paths. Stochastic differential equations, however, are far 
from stable under perturbations in the iterated integrals: recall the well-known examples of 
McShane IfTTII which give uniform approximations to Brownian motion where the limiting 
differential equation exhibits bias in the form of additional drift terms. (The explanation 
is that these approximations do not correctly approximate the iterated integrals of Brown- 
ian motion known as Levy's stochastic area.) At the risk of confusing the reader, even if 
is guaranteed that a sample path and its stochastic area are uniformly correctly approxi- 
mated, the limiting differential equation may still exhibit additional drift termsQ The point 
is that topology matters: uniform convergence needs to be replaced by a stronger notion 
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... now involving two or more iterated Lie brackets of the diffusion vector fields; see (4j for more on such 
subtleties. 
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of Holder (or ^-variation) rough path topology in order to use the stability results of rough 
path theory. 

Our key idea is to view the iterations of cubature steps, (Lyons and Victoir |9, The- 
orem 3.3] for instance), via an underlying random walk of the driving signal, Brownian 
motion plus Levy's area. The iterated cubature scheme corresponds precisely to stochas- 
tic differential equations in which the driving Brownian motion is replaced by k properly 
rescaled concatenations of the w, (say, chosen independently with probability Aj at each 
step). Thanks to the smoothness of the w,-, such a path has canonically defined iterated 
integrals; the "only" thing left to do is to establish weak convergence of this random walk 
to Brownian motion and Levy's stochastic area, in the correct rough path topology. It 
is then an immediate consequence of the continuity of the Ito map in rough path sense 
(i.e. as a deterministic function of path and area in rough path topology) to see that this 
entails the desired weak convergence result for path dependent functionals of the type 
f(S,:0<t<T). 

Weak convergence questions of this type were first discussed in E. Breuillard, P. Friz, 
M. Huesmann |2|. Unfortunately, the "Rough path Donsker" theorem obtained therein 
does not lend itself immediately to the present applications: a moment of reflection reveals 
that it would cover cubature with (1) equidistant steps and (2) in which the w,- are straight 
lines (Wong-Zakai!). Our strategy is thus to develop refined arguments that allow to cover 
the generic cubature setting as well as its recent variations (like Ninomya- Victoir). This 
leads, en passant, to a more flexible version of the Donsker theorem for Brownian motion 
on Lie groups in topologies considerably finer than the uniform one. 

The mathematical content - weak convergence of discrete structures to (Stratonovich) 
SDE solutions - should also be compared to the (typically Ito) diffusion limits of Markov 
chains (cf. Stroock and Varadhan [16. Section 11.2]), although we shall not pursue this 
point further here. 

The current paper uses many ideas and results from rough path theory, see Lyons [ 10 1 
and Friz and Victoir [5], which we primarily use for reference in this paper. For cubature 
on Wiener space, the authoritative reference remains Lyons and Victoir (5). 



2. Cubature on Wiener space and the associated random walks 

Let B - (Bf)re[o,i] denote a standard <i-dimensional Brownian motion on (Q, T, CFt)te[0,ih P) 
and B. = 5 l 2(-B)o, ) i- e > B is the Brownian motion enhanced by its Levy area. The geomet- 
rical setting of B is the Lie group G 2 (R d ), which can be defined as follows: let e\, . . . , ed 
denote the canonical basis of R d . Thene,®ej, 1 < i,j < d, forms a basis for the tensor prod- 
uct R d ®R d . Consider the algebra T 2 (R d ) := ReR< / eR<'®R fl ', which is understood as a step-2 
nilpotent non-commutative algebra, i.e., for Xi = z\ + x\ + a\,X2 — Zi + X2 + a 2 e T 2 (R d ) 
the product is given by 

xi ® x 2 = z\Zi + (zzxi + z\x 2 ) + (zia x + z\a 2 + X\ ® x%) . 

Consider g 2 (R rf ) c T 2 (R d ), the Lie-algebra generated by e\,. . ,,ed together with [e,,e ; ], 
1 < i < j < d, with the commutator defined by [x, y] := x ® y - y ® x, x, y e T 2 (R d ). The 
exponential map exp : T 2 (R d ) — > T 2 (R d ) defined by 



exp(x) := 1 + V — x® k = 1 + x + -x ® x 
£-l k\ 2 



ui k 



in the step-2 nilpotent setting, maps g 2 (R' / ) in a bijective way to the Lie group G 2 (R d ) := 
exp(g 2 (R rf )) c T 2 (R d ). 

This Lie group is highly relevant for rough path analysis, since it is the geometric setting 
of the enhanced Brownian motion mentioned before. Indeed, the T 2 (R rf )-valued process B 
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d a ' pt 

^B i t e i +'^ j I B i s odBie i ®e j =:S 2 (B)o,t, < f < 1, 



defined by 

■/ 

(1) B, := 1 + 

i=\ i,j=l 

lives in the Lie group G 2 (R d ), i.e., p(b, e G 2 (R d ), t e [0, 1]) = 1. In a similar way, we 
will consider the (step-m truncated) signature, see Friz and Victoir [5 1, 

m p 

(2) S m (B) y.= l + J) V odB i ^--odB i le k ®---®e k , 

k=li l ,...,i t e{l,...,d\ Jo ^-^' 

which takes values in the step-m nilpotent Lie-group G m (R d ) defined analogously to G 2 (R rf )o 
Consider the stochastic differential equation (in Stratonovich form) 

d 

(3) dX, = V (X,)dt + V Vi(X,) o dB\, 

i=i 

X = xq e R N . Here, V , V u . . . , V d : R N -» R N is a collection of smooth vector fields. 
A cubature formula on Wiener space is a random variable W taking values in the space 
C var ([0, l],R rf ) of continuous paths of bounded variation with values in R d such that we 
have 



| odB l ,\ ■ ■ ■ o dB\ = E I dWl 

J0<ti<-<t t <\ [Jo<t,<-<r t <\ 



■dW) k k 
-<t t <\ 



(4) 

■<fi<i 

Equation <(4j is supposed to hold for all multi-indices / = (i\,...,ik) e { 1 , . . . , d} k with 
k < m and all 1 < k < m, where m is a fixed positive integer, the order of the cubature 
formula. Moreover, we note that the paths of the process W are of bounded variation, 
therefore the integrals on the right hand side of (|4} can be classically defined in a pathwise 
sense. Notice that we do not use cross-integrals between time dt and the Brownian motion 
dB t . Therefore, a cubature formula in this sense can only be used to approximate SDEs 
with drift Vq = 0. We will cover the general case later in Section [3] 

Rephrased in terms of the (truncated) signature, equation (0| means that 

(5) E [S m (B) ,i] = E [S m (W)o,i] , 

where the expectation takes values in the algebra T m (R d ). Obviously, any cubature formula 
on Wiener space can be rescaled to a cubature formula on the interval [0, At], At > 0, by 
replacing W with the bounded variation path 

(6) d<jA,(W) : [0, Af] -» R d , sh \[AlW(s/At). 

On the level of signatures, this corresponds to applying the dilatation operator S^ : 

G m (R d ) -> G m (R d ), i.e., 

M 5 VA7(H,A, = <J VA# m W0,l). 

Remark 2.1. Note that the symbol 5 ^ has different meanings on both sides of the equa- 
tion: on the left hand side, it is a function from C( [0, 1 ] , R d ) to C( [0, At] , R d ), whereas on 
the right hand side it is the restriction to G m (R d ) of a linear map defined on the algebra 

T m (R d ) by 

6 VA7( e '. ® ' • • ® e k) '■= Afk/2e h ® ' • • ® e i- Q<k<m, i u ..., i k e{l,...,d}. 



Note that 5,„(B)o,. = S„,(B)o,., the Lyons lift of the enhanced Brownian motion, reflecting the fact that S m (B) 
depends uniquely and continuously on B - whereas B itself is not uniquely and certainly not continuously given 
by B. For instance, we could have chosen the Ito-integral instead of the Stratonovich integral. 
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Given a mesh D = {0 = to < t\ < ■ ■ • < t n = 1}, set Aft := ft - tk-i, k = \,...,n, and 
|D| = max* Aft. Moreover, let W(i), . . ., W(„) be independent copies of the cubature formula 
W. We define a random variable W fl : [0, 1] — > R ' taking values in the space of continuous 
paths of bounded variation by concatenation of the paths S^:(W(jc)) '■ [0, Aft] — > R d , 
k - \,...,n. Again, by well known properties of the signature (the Chen theorem, see for 
instance [2, Theorem 7. 1 1, Exercise 7.14]), this translates to the relation 

s m (w D ) oj = s m (*W w (i))) A;i ® • • • ® s m (s^r n ( w M)) 0M 

= S^(SJW (l) )o A )®---®6^(SJW {n) )o,i), 

where ® denotes the multiplication in the Lie group G m (R ). Finally, let X D denote the 
(pathwise ODE) solution of the equation 

(7) dX? = f j V i {XP)dW^, 

X® - xq. For a given function / : R^ — > R of interest, the method of cubature on Wiener 
space now consists in the approximation 

(8) E[f{X,)] = E [f{Xf)\ + O (|D| (m - 1)/2 ) , 

provided that certain regularity assumptions are satisfied, see (9), Ifl3l and Q. In partic- 
ular, the method provides an efficient numerical scheme, if W has been chosen in such a 
way that integration of (0 is "substantially simpler" then integration of the original Q, 
see QJ. If/ is smooth, (O holds even for uniform meshes. If / only is Lipschitz, how- 
ever, then Kusuoka [7 Theorem 4] shows that (JHJ holds provided that one takes certain 
non-homogeneous meshes. The goal of this paper regarding cubature is to show that con- 
vergence even holds for (reasonable) functionals / depending on the whole path (X t )o< t <i. 

Example 2.2. The cubature formulas in [9| are discrete random variables W taking values 
in the space of continuous paths of bounded variations. That is, fix k paths of bounded 
variation o>i, . . . , a>k '■ [0, 1] — > M. d and positive real numbers A\, . . . , A^ with A\ +■ ■ ■ + A/ i — 
1. Then W is the random variable taking values in {<wi, . . . , a>k) with P(W — coj) — Aj, 
j — 1 , . . . , k. In all the concrete cubature formulas constructed in [9 1, the paths tOj(-) are, 
in fact, piecewise linear. 

Example 2.3. We can even interpret the Wong-Zakai approximation as a cubature formula 
on Wiener space (of order m — 3; the resulting convergence in © has then weak order 
(m — l)/2 = 1, precisely as the usual Euler scheme for Ito differential equations). Indeed, 
choose W as the linear path W t — tB\. We note that W can be realized (for any mesh 
D) by choosing 6 ^(W^Xs) = s(B tt - B, k _ x ), because W (k) (t) := -^= k (B, t - B,^) has the 
same law as W and all the W(k) are independent. Concatenation of these paths precisely 
gives the piecewise linear approximation of B with nodes in D. 

Example 2.4. Ninomiya and Victoir [ 13 1 construct a cubature formula of order m — 5 in 
the following way. Let Abe a Bernoulli random variable (taking values ±1 with probability 
1 /2 each) and let Z 1 , . . . , Z be independent standard normal random variables. Set s — 
l/(d + 1). For a> € Q, W(co) is defined by the following formula. If h.(a>) — — 1, we define 
W(oj) to be the piecewise linear path with 



W\cjXs) 



\/e, se [0,s/2], / = 0, 

Z'(w)/e, s e]e/2 + (i - l)e, e/2 + re], ie{l,...,d], 

1/e, se]l-6/2,l],j = 0, 

0, else. 
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If A(oj) — 1, W(co) is similarly defined by 



W'( W )(5) = 



1/e, s£[0,e/2],! = 0, 

Z'(oS)/s, s e]e/2 + (d- i)s, e/2 + (d-i+ l)s], i e {1, . 

1/e, je]l-e/2,l],( = 0, 

0, e/se. 



..,<*}, 



77z/s means, we subdivide the interval [0, 1] into d + 2 subintervals 



K 


u 


s 3e 

2'T 


u- 


■ u 



I_3f,l_f 

2 2 



2 



On eac/i of these subintervals, W(co) is constant in all components albeit one, which is 
linear. In particular, W(aS) is again piecewise linear. 

At this stage, we would like to remark that we could replace the Gaussian random 
variables Z' by discrete random variables having the same moments of order up to five. 
Then we would obtain a special case of Example \2.2\ - albeit for the non-standard choice 
ofW°, see Section\3\below. 

Let us now turn our attention to Donsker type results: for a fixed sequence of meshes 
T> n with \T> n \ — > we wish to study the corresponding sequence of paths in G m (K. d ), i.e., 
we study 



M<n,4-(<n,L, 11 



By a Donsker theorem in rough path topology for the sequence of cubature formulas W v " 
we understand the statement that 



(9) 



S m (w D ") >S m (B) , 



weakly with respect to or-Holder rough path topologjQ for some a e (1/3, 1/2) andm > 2. 
(In fact, elementary results of rough path theory imply then that it suffices to consider 
m — 2. Also, the claimed convergence will actually be established for all a < 1/2). As 
a justification for calling the convergence stated in ((9} a Donsker theorem, consider the 
following random walk. Let us again fix the mesh T>„ = {0 = to < ■ ■ ■ < t n = 1}. (We 
only take n as the number of sub-intervals for the grid D„ for more convenient notation. 
The mathematics would, of course, work in precisely the same way, if the size of D n was 
completely arbitrary, as long as \D n \ — * 0.) Define 



(10) 



$ = s m(s ^rSW {k) )) Qhit = 5^(S m (w (k) X\), 



a random variable taking values in G m (R, d ). Note that £jj? = 6 jA7"(£ot))> wriere %(k) is an 
independent copy of 5 m (W)o,i. Next define the G'"(R. fl ')-valued, finite random walk S?, 
k = 0, . . . , n, by Z" = 1 and 

"*+l _ "* w $k+l> 
where 1 is the neutral element of G m (R. d ). Since 



s 4 w °X, t -^ fe =°- 



5 m (W B " J is, indeed, a path in G m (R d ) obtained from the random walk S" by (possibly 
random) interpolation. This gives the link to the classical Donsker theorem as well as to 
the paper of Breuillard, Friz and Huesmann [2|. Let us rephrase their Theorem 3 for the 
current setting. 



I.e., the ff-Holder topology for functions taking values in the metric space (G 



, see (5] Definition 



5.1,9.15] 
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Proposition 2.5. Let W be a cubature formula on Wiener space of order m — 2 with finite 
moments of all orders in the sense that 

V?>l:£[||S 2 (W)o,i|| ? ]<c», 

where \\-\\ denotes the Carnot-Caratheodory norm on G 2 (R. d ), see below. Moreover, assume 
that W is chosen in such a way that for every a>, S 2( W(lj))o,- is a geodesic connecting 1 and 
S 2(W(ui))o t i. Choose uniform meshes T>„ — I - k = 0, . . . ,n\. Then the Donsker theorem 

holds in rough path topology, i.e., S 2 (w D ") converges to B in C°*- m ([0, 1], G 2 (R d )), 
for every a < 1/2. 



Recall that the Carnot-Caratheodory norm is defined by 



infj f \dy\ y e C 1_vflr ([0, l],R d ), 



S 2(7)0, 



The infimum is actually always attained, and can be parametrized as Lipschitz-continuous 
path with constant speed, i.e., S2(y)o,r = f ||5'2(y)o,i|| for < t < 1 and a minimizing 
path y, see {5\ Theorem 7.33]. As a homogeneous norm, the Carnot-Caratheodory norm is 
equivalent to the simpler norm 

||x|| 2 := max(|jc| , |a| 1/2 ) , x = 1 + x + a e G 2 (R d ), 

see J3J Theorem 7.45]. 

Remark 2.6. In [2 Theorem 1], the moment condition is relaxed, which gives weak con- 
vergence in a-Holder norm for all a < a*, for some a* < 1 /2, which is related to the 
relaxed moment condition. In this paper, we shall always assume existence of all the mo- 
ments. We note, however, that we could also relax this assumption, obtaining a similar 
result. 

3. The main result 

Usually, a cubature formula W will not satisfy the conditions of Proposition ^. 51 even if 
we only choose uniform meshes, because the corresponding interpolation S \W D " ] of the 
random walk H" will not be geodesic. Moreover, if we want to treat functions / which are 
not smooth, then we have to choose non-uniform meshes with h = ^ for some y > m — 1, 
see 1 7 1. Therefore, we want to generalize Proposition 12. 5 l in two directions. We want to get 
rid of the condition of geodesic interpolation, and we want to generalize to non-uniform 
meshes. Fortunately, the first generalization is simple, at least for the cubature formulas 
actually suggested in the literature, see Example 12.21 12.41 and also for the Wong-Zakai 
approximation given in Example 12. 3 1 The second generalization, however, requires us to 
change the method of proof as compared to |]2]. 

It is natural to impose some restriction on the behavior of S \W D " ) between two nodes 
of the random walk. Indeed, we have to rule out "loops" which approach infinity. 

Assumption 3.1. The cubature formula W takes values in the Cameron-Martin space "H 
(of paths started at 0) and the Cameron-Martin norm has finite moments of all orders, i.e., 
for every ieN 

,1 \*/2l 



E[|IW|iy = £ U \W(s)\ 2 ds 



This assumption is both natural (a general continuous path of finite 1 -variation would 
not be in the (l/2-e)-H61der support of the Wiener measure!) and satisfied by all (piecewise 
linear!) cubature formulas used in practice. We look at this in some detail in 
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Example 3.2. Assume that the cubature formula W is piecewise linear, i.e., there is a posi- 
tive integer I and there are d-dimensional random variables F\,...,Fi with finite moments 
of all orders and a mesh = sq < ■ ■ ■ < S{ — 1 such that 

W s = Fu si-i < s < si, 1 < / < (. 

This immediately implies Assumption ^. 71 

Our main theorem is (for conclusions to cubature see Corollarv l3.5l below): 

Theorem 3.3. Given a cubature formula W of order m > 2 such that W\ and the cor- 
responding area A\ have finite moments of all orders and Assumption l3.il is satisfied. 
Then Donsker's theorem holds in rough path topology for any sequence D n of meshes 
with \T> n \ -> 0, i.e., 

Siiw®") >S 2 (B) , 

in C°' a - H "'([0, l],G 2 (R d )), for every a < 1/2. 



The natural conclusion from Theorem 13. 3 1 would be a weak convergence result for the 
cubature-approximation of the SDE d3} to its true solution on path-space. A little care 
is necessary, however, because we have ignored the drift Vo in the SDE, i.e., our driving 
signal is a pure Brownian motion and does not include time. The classical approach is to 
add another component to both the Brownian motion and the approximating cubature paths 
by setting B® := t, Wf := t and then require the moment matching condition d4} to hold for 
all iterated integrals, where the multi-index (i\, . . . , i%) now varies over {0, 1, ... , d} k , i.e., 
where we also consider mixed iterated integrals of Brownian motion and time t. Due to 
the scaling of Brownian motion "dB t a yfdt", it is only necessary to impose the moment 
matching condition for multi-indices (i\, . . ., i^) with k + #{j\ij - 0} < m to get weak 
convergence with rate 2yl . 

However, the Ninomiya-Victoir scheme does not fall into this class, because we have 
seen in Example l2.4l that they do not choose Wf = t. Therefore, we want to generalize the 
above considerations slightly. Let h : [0, 1] — > R be a deterministic, uniformly Lipschitz 
path with h(Q) = and h(\) - 1. This setting obviously includes the drift-component of 
the Ninomiya-Victoir scheme. We define the path h W by h W' t := W' t for i = \,...,d, and 
h W® := h(t). As usual, we set B® := t. We assume the usual moment matching condition to 
hold, i.e., 

E[s m ( h W) ,i]=E[S m (B\i], 
where S m ( h W) is the step-m signature of the path h W, more precisely 

m p 

M' ,W )oi=X Z d h W? i --d h WZe k ®--®e k . 

*=U (j 1 ,...,! t )€{0,l,...,d) 
k+#{j:ij=0}<m 

Analogously, the signature of the Brownian motion above is understood as the signature 
of the now R^'-valued process B. (This notation is ambiguous. In the following, the 
symbol B will usually denote the Revalued Brownian motion. We only mean the extended 
R° f+1 -valued process if specifically indicated.) We note that the signatures of the (d + 1)- 
dimensional processes take their values in a stratified Lie group denoted by G'"(M. d ). In a 
similar fashion as above, we obtain - by rescaling and concatenation - a stochastic process 
h W along a grid T>. Of course, we have to use a different rescaling for the component 
h W°. Indeed, following the construction in Section|2] we define S^( h W) : [0, At] -> R d+1 
by Stfjd'Wy, = VAfW; /A , for i = 1,.. .,d, as before, but S^( h W)° s = Ath(s/At). We 
continue to construct h W v by concatenation. 
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By rough path theory (continuity of Young pairing, e.g. ||5] Section 9.4.4]) we also 
obtain weak convergence in path-space for the extended process h W to the extended, Re- 
valued Brownian motion, which even holds for any truncated signature - not only for the 
step-2 signature. 

Corollary 3.4. Let W be cubature formula on Wiener space of order m > 2 satisfying 
Assumption \3. 1 1 and such that W\ and the corresponding area A\ have finite moments of 
all orders. Then Donsker's theorem holds in rough path topology for any sequence D„ of 
meshes with \D„\ — > 0, i.e., 

Sn(''W d ") >S N (B)o, 

in C°' a - m ([0, 1], Gf (R d )),for every a < 1/2 and any N > 1. 
Moreover, we define h X D as the solution to the (random) ODE 

(11) d h Xf = V (*S?) dh v " (t) + J] Vi ("Xf ) dWp'\ 

where h D " := h W v, "°. Then we have weak convergence of h X :D " to X on path-space. 

Corollary 3.5. Given a bounded, continuous functional f : C°' a ~ Hol ([0, 1],!^) — > R, and 
assume that W, h and D„ satisfy the assumptions of Corollarv \3.4\ Then we have 

E\f( h X D ")] >E[f(X)l 

where X denotes the path (X,), e [o j i] of the true solution of the SDE Q and h X D " denotes 
the pathwise solution of the ODE ( llll l. 

Proof. We interpret <JVj as a rough differential equation, i.e., for a given (rough) path a> e 
C°' a - m \[0, 1 ] ; G^W 1 )) with 0) = n\ (to) we define n(<o), := y, by 

d 

dy, = V Q (y,)dJ! + J] V t (y t )dcJ t . 

i=i 

In particular, we have h X°" = n(s 2 ( h W D ")o, .) andX = 7r(5 2 (B)o,)- By (1 Theorem 10.26], 
the map a) h> 7r(w). is a continuous map from C aff " H81 ([0, 1]; G\(R d )) to C 0>Qr - H61 ([0, 1]; R w ). 
Thus, CorollaryElimplies weak convergence of ''X " = tt (5 2 ( /, W°")o,) to X = 7r(5 2 (fi)o,) 
in C°^ H81 ([0, 1];^). a 

Remark 3.6. 5/nce f/ie a-Holder topology is stronger than the usual uniform topology 
given by the supremum norm, Corollarv \3.5\ in particular holds for all bounded functionals 
f which are continuous in the uniform topology on path space. In the case of unbounded 
continuous functionals, convergence can still be guaranteed provided that some uniform 
integrability property holds. (Of course, in the case of call-option type derivatives, one 
could also try a relevant put-call-parity.) Finally, in the case of barrier options, the payoff 
functional is often continuous apart from a set of measure zero on path space. Naturally, 
non-continuities on null-sets do not hinder weak convergence of the cubature method. 

4. Random walks with independent, non-identically distributed increments 

In this section, we prepare the main ingredients of a proof of Donsker's theorem for 
random walks with independent, but not identically distributed increments on the Lie group 
G := G 2 (R rf ). More precisely, let £ be a random variable with values in G with finite 
moments of all orders. We shall denote the components of £ in the basis of g := g 2 (R"') 
given by e,, 1 < i < d, together with [e,, ej], 1 < i < j < d, by X' and A'-i, respectively, i.e., 

' d \ 

£ = exp ^X< e , + JV'-' / [ e; , e j] 

V i'=l i<j 



CUBATURE ON WIENER SPACE: PATHWISE CONVERGENCE 9 

Thus, the condition that £ has finite moments of all orders simply means that all the real 
random variables X', AS have finite moments of all orders q > 1. Moreover, we assume 
that £ is centered, i.e., 

£[X'] = 0, i=l,...,d. 
Let us fix a mesh T>„ = {0 = to < ■ • ■ < t n = 1}. For n independent copies £(i), . . . ,£(„) 
of f, define ^ = S ^-(if^)) and the corresponding random walk 

For use in the next lemma, let us define the coordinate mappings x' (mapping ieGto the 
component of log(^) with respect to the basis element e,) and jc' j (mapping x e G to the 
component of log(.x) with respect to the basis element [e,-, e 7 ]), 1 < i < d, i < j < d. As 
usual, the corresponding vector-fields (i.e., basis of the tangent space) are denoted by -^ 
and -Aj, respectively. 

Lemma 4.1. The above random walk satisfies the central limit theorem, i.e., Ej| converges 
weakly to the Gaussian measure with infinitesimal generator 

y a uJ_ + i_y bU dd 

Zj dx'-J 2^ dx' dxJ 

KJ HJ 

where a" := E[A iJ ] and b iJ := Cov(X',X J ')fl 

Proof. The result is well-known in probability theory on Lie groups, see, e.g., 1141 . We 
verify that the system of probability measures fi,,^ - (gT)*P, i.e., [i n ^ is the law of £?, 
satisfies the conditions given in [14, Theorem 3.2], namely: 
(i) sup„ Y1U\ Ic \\x\\ 2 Vn,k(dx) < oo; 
(ii) fi ni k is centered in the above sense; 
(iii) for every 1 < i < j < d, there is a number a'-' e R such that 



" r 

= ]mt r. 

k=\ U(j 



x l ' J (x)n„jc(dx); 



(iv) for every 1 < i, j < d, there is a number b'-> e R such that 
b' J ' = lim y I x i (x)x^(x)ii n i i {dx); 



*=i 



(v) lim^ Y!U Jf Wfee IMfWW*) = for all e > 0. 
By homogeneity of the Carnot-Caratheodory norm, we have 



Thus, 



J \\x\\ 2 M „, k (dx) = E [||<5v5r^)|| 2 ] = At k E [ll^ll 2 ] • 

n p n 

£ I ||x||Vm(^) = J] Af,£ [ll^l 2 ] = £ [ll^l 2 ] , 



t=l ^° jfc=l 

and the supremum over n is obviously finite, settling (i). 

(ii) is satisfied by assumption on £. Regarding (iii), note that x' J (^p = x'^{5^^{^)) - 
AtkA l, i, where equality is understood as equality in law. Therefore, (iii) is satisfied with 
a l 'J - E{A''i) < oo. A similar argument shows that (iv) holds with V-' - Cov(X' ,-X 7 ). 

For the proof of (v), we again use homogeneity of the Carnot-Caratheodory norm. In- 
deed, we have 

f IMrWdx) = & k E [i ]6 , M[ ( VM^II) \\f\\ 2 ] , 



The statement means that there is a semi-group (h,),>q of probability measures on G having the above 
infinitesimal generator and /i i is the limiting distribution of EJj . Moreover, this semi-group is Gaussian in the 
sense that lim^o [Mt(P \ U) = for every neighborhood U of the neutral element of the group G. 
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implying that 






||*|| n n ,k(dx) < E 



v.')„r 




V^M' 



by the Cauchy-Schwarz inequality. Now, the right hand side converges to zero by integra- 
bility of ||£|| and |D„| -> 0, for every fixed e > 0. n 

Remark 4.2. If g is the step-2 signature of a cubature formula of degree m > 2, then 
a iJ = E[AU] = 0, 1 < i < j < d, and, moreover, V'> = Q.qm{X 1 ,X') = 6 U . Thus, the 



generator of the limiting Gaussian measure in Lemma \4~l\ coincides with the generator of 
the Brownian motion on G, i.e., with the generator o/B. 

Next we state a moment estimate, which will enable us to prove tightness of the family 
of interpolated random walks in rough path topology. 

Proposition 4.3. For every p e N, p > 1 we can find a constant C independent ofk and n 
such that 



| s «||4p 



<Ct 



2/> 



Proof The proof heavily relies on Burkholder's inequality, see O. Recall that the discrete 
time Burkholder inequality establishes the existence of constants c p , C p for 1 < p < oo such 
that for every p-integrable real martingale Y„ and any n e N we have 



c p sup \\S„\\ L „ < sup ||F„|| < C p sup ||5„| 



Li' 



where, setting Yq :- 0, S„ :- ^Y!l=\^k - l*-i) 2 is the square root of the quadratic vari- 
ation of Y. By choosing Y„+i = Y„ for / > 0, this immediately implies the corresponding 
finite version 



(12) 



c P lls„|| LP <||y„|| LP <cj|s„|| 



By equivalence of homogeneous norms, see, for instance, [5, Theorem 7.44], we can re- 
place the Camot-Caratheodory norm ||-|| on G 2 (R d ) by the homogeneous norm 



||x|| 2 := max(|^aog(x))| , A /|^ 2 (log(x))|) < |^(log(x))| + ^\7T 2 (log(x))\, x e G 2 (R d ), 
where n\ and n2 denote the projection to the first and second level components of x, i.e., 



when x= I + x + a e G (R ), then ^i(x) - x e W and 7r 2 (x) = a e 



Thus, the 



assertion of the proposition is equivalent to the existence of a constant C (only depending 
on p) such that 



(13) 
(14) 



Fi 



oog^))! 4 ' 

^(log(Ep)| 2 " 



< Ct 



i P 



<Ct 



2/' 



We start by proving (Q~3). By the Campbell-Baker-Hausdorff formula, we have 



Fi 



(log(H£))| 



d ( k \ 

z ^ « 


d 

/=i 


k 



Here, X\ - y/AfiXL and C is a constant, which does neither depend on the partition D n 
nor on k. For the remainder of the proof, we will use this symbol for constants that may 
vary from line to line, but do not depend on £)„ or on k. This implies that 

4pi 



^(log^))! 4 " <C^E 



I*; 

/=! 
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Now we apply Burkholder's inequality to the martingale Y^ - Y!i=i X\ with the exponent 
4/7 to get 



|7n(log(H£))| 4p ]<C^£ 



1=1 



2/' 



ctf2 E 



i=i 






/=1 



f* 



2/' 



Noting that the sum inside the expectation is a convex combination, we apply Jensen's 
inequality for the convex function x 2p and get 



d k 



;=i 1=1 



, l( iog (E «))|^ < c? Y,Y,t e [<V1 = c Z £ K>>*] 



z 



#> 



which is of the form required in ([T3V 

For (fl4l . we again start with the Campbell-Baker-Hausdorff formula and get 



h(log(H£))| 



■="i\\l 2 ^ 



= £ 



Z Z^4 Z (KK-KK) 



1<i<j<dll=1 



1<l,<h<k 



(15) 



<C £ £ 

1 <i<j<d 
l</'</'<d 



Z4'M Z (44-44) 



2/' 



2/' 



/=1 
* 



1<l,<h<k 



Z^ 

;=i 



2/' 



z 



xj xj 

h h 



1<Ii<h<k 



2p 



+ E 



Z 44 



l<l,<h<k 



2p 



Now fix some 1 < i < j < d. Again by Jensen's inequality, we have 



(16) E 



z^ 

/=1 



2/' 



?* 



V — A u 

6 1 f * (0 



2/' 



< r 



2/' 



A/v 



7T £ 



/=! 



\A' J 
\ A (l) 



2/' 






flp 



?■ 



On the other hand, note that 



(h-1 \ 



Z 44=Z4Z4 



l<Zi<fe<* 



/ 2 =l V/i=l / 



is a martingale (indexed by k). Thus, Burkholder's inequality (TTZt for the exponent 2/? 

gives 





2p" 




k 


ffe-i V 


p 




Z 44 




<CE 


Z^) 2 


Z4 




= Ct" k E 


/l</ 2 St 






'2 = 1 


Wi=i / 







Z^<<,' ! [Z4 



Z4 

Ui=l / 



Now we again apply Jensen's inequality and then the Cauchy-Schwarz inequality, and 
obtain 



z 



X) X] 

h h 



1<h<h<k 



2p 



k 



n V "tin 

^Zt* 



k 



h 



(4)) 2 



(h-\ \ 
Z4 

V/, = i ) 



(17) 



* <*£*£(* 



ii=i 



\ A (h) 



4p-l\l/2 ' 



k-1 

Z4 

ii-i 



4pi\ 1/2 
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By applying Burkholder's and Jensen's inequalities for a final time, we get for the left-most 
term in the above inequality 






Ji-i 



M< 



< Ct 2 "E 



fe-l A 

V^ Af/, i ■) 


2p- 






<Ct] 


■P 

2-1 


^ Aft, 
> — -£ 


Ir l 4;/ 



< Ct] p _ x E 



\y 



|4p 



(Di 



Inserting the last inequality into ( TTTl i. we obtain 



Z n4 



l</l</ 2 <£ 



2,, 



^g^t,(+;,,P])>N1f 



t 2p 

Together with ([T3T > and (TT6b this shows ( TBI , and the proposition follows. 



5. Proof OF THE MAIN RESULTS 

Analogously to [2, Theorem 1] we can now state our 

Theorem 5.1. Let D n — {0 — to < ■ ■ ■ < t n — 1} be a sequence of meshes with \D„\ — * 
and let E" = (EjM be a centered random walk in G 2 (R, d ) along the mesh (as defined 
in Section |3, whose increments have moments of all orders. Additionally, we impose 
E{7tz(ffu) — 0. Define a sequence E of stochastic processes with values in G (R ) fey 

— n 

E r = "Bl^for k = 0, . . . ,n and by geodesic interpolation for t e [t^, tk+i]. Then 



^B 



in C a '°- Hel ([0, 1], G 2 (R d ))for all a < 1/2. 



Proof Kolmogorov's tightness criterion, see, for instance, lUD . implies that S is tight (in 
C°'^ H81 ) provided that for any u, v 6 [0, 1] 



(18) 



[/ — n — n \ a ~\ 
d \-v^u) \<C\V-U\ 



where a, b, c are positive constants with a < b/a and d denotes the Carnot-Caratheodory 
distance defined by d(x,y) - \\x~ 1 y\\. We choose a — 4-p and b = 2p — 1, then Proposi- 
tion l4.3l implies that ( [T"8l holds for u,ve D„. For arbitrary v < u, assume that t ; < v < f,+i 
and tj < u < tj+\ (for f,-, f,+i, fy, f/+i e £)„). Using (by the geodesic interpolation) 



(19) 



'PC*.) 



d K • ^ ) ' d K> *u ) = -TT-d (^ , ^, ) 



Af,- H 



Af/+i 



and the triangle inequality, we obtain 

[I — n — n\a~\ 
d^ u )]<c 



fc+i - V 



Aft+i 
<c\u- v\ 2p , 



(f, +I -y) 2 " + (f i --f, +1 ) 2 ^ + 



Af 



7+1 



IP 



(u-tjfP 



for some constant c only depending on p. 

This shows that the sequence of stochastic processes S is tight in rough-path-topology. 
Moreover, Lemma |4~T1 shows that the finite-dimensional marginal distributions of S con- 
verge to those of B. Thus, we obtain the theorem for a < j-j^ and, with p — > oo, for any 

a < 1/2. □ 
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Proof of Theorem \3.3\ In Theorem l5. ll above. we have already proved our main result for 
N - 2 and for the special case of geodesic interpolation, i.e., for the case that S 2 (W D,, J , 

t e [0, 1], provides a geodesic interpolation between the grid points, i. e. , between^ (W fl ") 
tk € T)„. We note that the extension of the result to Af > 2 is immediate since the Lyons 
lift S n : G 2 (R d ) — > G N (R d ) is continuous in rough path topology, see, for instance, [5, 
Corollary 9. 11]. 

We need to give a proof for the case N - 2 but without geodesic interpolation. As 
before, assume that we are given t; < v < f;+i, tj<u< tj+\, with f,-, tj, f !+ i, tj+i e T> n . For 
any path a> e < H, the Cameron-Martin space, we have, see ||5] Proposition 15.7], 



where I 



I l-var;[s 



Ml-var;[v] ^ Vl f - s \ IMI«;[ S ,fl , 

,j denotes the first variation of a path restricted to [s, t] and ||-||^ ; [ S ,] denotes 



the Cameron-Martin norm, likewise restricted to [s, t], i.e., 



Notice that 



\W 



D„ I 



IMI«;[V] = J l^ M 

[ '-" ] J, I I A^i J,, 



' du. 



./'+' 



rff < 



rV 1 

w,' 

Jo 



IW;[f 

Therefore, we can bound 

1 1 l-var;[ry,w] 

< ,/u=Tj\\W\\f<- 



dt 



By Assumption l3.ll we get 



(20) E 
and similarly 
(21) 



d ( s * k% "^ K°)o,«n - (m _ ° )2 " £ [ |iw "«] - c(m _ f;)2 "' 



^Na^iaf 



< C(t i+l - v) 



2/; 



Therefore, we can show tightness of the sequence of processes S 2 (W 23 "} as in the proof 
of Theorem l5. ll above. □ 

Now we finally return to the Donsker theorem for cubature paths with an adjourned, 
Lipschitz component h. 

Proof of Corollar\ \3~4\ Let h D " := ''W '" denote the 0-component of h W D ". Moreover, 
let id : [0, 1] — > [0, 1] denote the identity, id(f) = t. We note that h " converges to id in 
C°^ _H81 ([0, 1];R) for any /? < 1. Indeed, let < t < 1 and let 1 1 < t < t M denote the grid 
points closest to t. Then apparently 



hr n {f) = U + Atj+ih 



t-tj 

At, 



/+! 



For h(0) = 0, we get \h{t)\ < L\t\, where L is the Lipschitz constant of h. From this one can 
easily conclude that yfr " - idj^ < (1 + L) \D n \ and h D " converges to id uniformly on [0, 1] 
with uniform Lipschitz bounds for \D„\ — > 0. This implies the convergence in^S-Holder 
topology for any (5 < 1 . 

By this result together with Theorem 13. 3 1 for the convergence of S 2 {W), we can im- 
mediately conclude that S 2 { h W D ") converges to S 2 (B)o, in C°'^ Hel ([0, 1]; G?(R d )). By 
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continuity of the Lyons lift S2W l_> Sn(x) in rough path topology, the statement of the 
corollary follows. □ 
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